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Abstract 

We present a systematic perturbative construction of the most general metric operator 
(and positive-definite inner product) for quasi-Hermitian Hamiltonians of the standard 
form, H = \p 2 + v(x), in one dimension. We show that this problem is equivalent to 
solving an infinite system of iteratively decoupled hyperbolic partial differential equations 
in (1 + l)-dimensions. For the case that v(x) is purely imaginary, the latter have the 
form of a nonhomogeneous wave equation which admits an exact solution. We apply 
our general method to obtain the most general metric operator for the imaginary cubic 
potential, v(x) = iex^ . This reveals an infinite class of previously unknown CVT- as well 
as non-CPT-inner products. We compute the physical observables of the corresponding 
unitary quantum system and determine the underlying classical system. Our results for 
the imaginary cubic potential show that, unlike the quantum system, the corresponding 
classical system is not sensitive to the choice of the metric operator. As another application 
of our method we give a complete characterization of the pseudo-Hermitian canonical 
quantization of a free particle moving in R that is consistent with the usual choice for the 
quantum Hamiltonian. Finally we discuss subtleties involved with higher dimensions and 
systems having a fixed length scale. 
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1 Introduction 



In Quantum Mechanics the unitarity of the dynamics dictates that the Hamiltonian operator 
must be Hermitian (self-adjoint). Therefore a non-Hermitian operator H cannot serve as the 
Hamiltonian for a unitary quantum system unless the Hilbert space is endowed with a different 
inner-product which renders H Hermitian. Under rather general conditions, one can show that 
a necessary and sufficient condition for the existence of such an inner product is the reality of 
the spectrum of H , pQ . The most general inner product that makes H Hermitian has the form 

<V>„ + :=<-b7+|->, (1) 

where (-|-) is the defining inner product of the original (reference) Hilbert space H, r] + : H — > H 
is a positive-definite operator satisfying the pseudo- Hermit icity condition [2] 

#t = V+ H V -\ (2) 

and W : 7i — > H is the adjoint of H. 1 

The operator rj + , which is sometimes referred to as the metric operator, is the basic ingre- 
dient of a quantum theory based on the Hamiltonian H. It determines the inner product (JJJ) 
of the physical Hilbert space 7i p h ys and specifies the observables O of the theory, |U 03 El ■ 
Furthermore, it provides the necessary means to obtain an equivalent description of the system 
using a Hermitian Hamiltonian h acting in 7i. This in turn allows for the determination of the 
underlying classical system, jHJ El • 

An important property of the metric operator, is that it is not unique [HI El HOI ITT] . For 
example, the CPT-inner product proposed in can be conveniently expressed in the form (0) 
for a particular choice of r] + , (HI US]. As such, the CPT-inner product is not the only allowed 
choice for the inner product of the physical Hilbert space; it is merely an example of an infinite 
set of possible choices. 

In jSJinilZl we have constructed specific metric operators for different PT-symmetric systems 
and examined the physical consequences and the classical limit of the quantum systems they 
define. In the present article we intend to address the issue of the non-uniqueness of the metric 
operator. In particular, we will offer a perturbative construction of the most general metric 
operator for the standard Hamiltonians in one dimension, and explore the consequences of 
adopting different metric operators for a given Hamiltonian operator. 

Our method which is based on a variation of the approach of jU] has the following two 
advantages, (i) It does not involve using a particular ansatz for each choice of the potential. 
Instead, it employs the well-established properties of (1 + l)-dimensional hyperbolic partial 
differential equations. This makes it applicable for all potentials, (ii) It is not restricted to 
producing the CPT-inner products; it yields the most general inner product that is capable of 
restoring the unitarity. In fact, as we will see for the particular example of the imaginary cubic 
potential, it reveals a new set of CPT-inner products that has escaped the analysis of [Tlj . 

1 For any densely defined H, is the unique operator satisfying ■) = (£[■ |-), 
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2 Perturbative Construction of the Metric Operator 

Consider a Hamiltonian of the form 



H = H + eH 1 , (3) 

where e is a real perturbation parameter and H and Hi are respectively Hermitian and anti- 
Hermitian e-independent operators. Suppose that for sufficiently small values of e the Hamil- 
tonian H has a real spectrum and a complete set of eigenvectors, so that a positive-definite 
metric operator r] + exists jj. Such an operator has a well-defined Hermitian logarithm. We 
shall let Q := — In 77+, or alternatively 

V+ = e~ Q , (4) 

and employ (J2J) to obtain a perturbative expansion for Q. 

The condition that for e = 0, the inner product (-, -) v+ of (JTJ) reduces to (-|-), i.e., as e — > 
we have 77+ — > 1, justifies the following perturbation expansion for Q. 

00 

Q = ^2Qj^, (5) 

3=1 

where Qj are e-independent Hermitian operators. A convenient tool for calculating the opera- 
tors Qj is the Baker-Campbell-Hausdorff formula, 

00 1 

e^He^=H + Y^^[H,Q] k , (6) 



k=l 



where [^Q], := [H,Q] and Q] fc+1 := [[#,Q] fc ,Q] for all k E Z+ := {1,2,3, •■•}. As we 
show in Appendix A, we can use (JHJ) to obtain the following infinite set of operator equations 
for Qj. 

[H ,Q j ]=R j , (7) 

where j G Z + is arbitrary and 

f -2/Ti for 7 = 1 

Rj := < ™ (8) 

{ Efe= 2 9* for 3 > 2 , 

m=l n=l \ / 

3y := £ [[•■■[^o,Q Sl ],g fl2 ],--- (10) 

si H h s/c = i 

It is not difficult to see that the anti-Hermitian operators i?j are uniquely determined by Hi, 
Qii Qii • •• ,Qj-i- This shows that (JJJ) may in principle be solved iteratively. For example 
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setting j = 1, 2, 3, 4 and 5, we respectively find 



[H ,Qi] = -2H 1} (11) 
[H Q ,Q 2 ] = 0, (12) 
[H ,Qs] = hH ,Qi] 3 =-l[H l ,Q l ] 2 , (13) 



12 u '^ x " 6 
1 

12 
1 



Wo,Qi] = ^([[Ho,Qi} 2 ,Q 2 ] + [[[H ,Q 1 },Q 2 },Q 1 } + [[H ,Q 2 },Q 1 } 2 ) 



~([[H 1 ,Q 1 ],Q 2 ] + [[H 1 ,Q 2 ],Q 1 ]), (14) 

D 

[H , Q 5 ) = — [Ho, Qi] B + ^ ([[Ho, Qi], Q 2 ] 2 + [[H , Q2], Qi] 2 + [[Ho, Qa] a , Qi]+ 
[[Ho, Qi] a , Q 3 ] + [[[Ho, Qi],Q 3 ],Qi] + [[Ho, Q 3 ], Qi) 2 ) 
= ^[H 1 ,Q 1 ] 4 -^([H 1 ,Q 2 ] 2 + [[H 1 ,Q 1 ],Q 3 ] + [[H 1 ,Q 3 ],Q 1 ]), (15) 

where we have made use of the fact that qt vanishes for even values of k. 

Comparing (fTT]) - ()15|) with operator equations obtained in jHj for Qj, we see that the 
choice of the CPT-inner product corresponds to taking Q 2 — Q± = ■ • • — Q 2 j = ■ ■ ■ = 0. 
Furthermore this choice restricts Q 2 j+i to the solutions of (JJJ) that are even in x and odd in p. 
We will refrain from imposing such restrictions and try to obtain the most general solution of 
0. 

Setting H = P 2 + Vo(x), where Vq is a real- valued function, adopting units in which 
m — H — 1, and evaluating the matrix-elements of both sides of (j7J) in the x-representation, we 
find 

[-dl + d 2 y + 2v (x) - 2v (y)](x\Q j \y) = 2(x\R j \y). (16) 

This is a nonhomogeneous hyperbolic partial differential equation in (1 + l)-dimensions. In 
particular, for the cases that v = 0, i.e., H = \p 2 , it reduces to the (1 + l)-dimensional 
nonhomogeneous wave equation, 

[-d 2 x + d 2 }(x\Q 1 \y) = 2(x\R 1 \y), (17) 

whose general closed-form solution is well-known [TH] . 

In summary, we have reduced the problem of determining Qj to the solution of a set of itera- 
tively decoupled and well-known partial differential equations. In fact, because these equations 
only differ in the form of their nonhomogeneity, we can express their general solution as 

(x\Q J \y) = f 3 (x,y) + Q 3 (x,y), (18) 

where fj is the general solution of the homogeneous equation 

[-81 + d 2 y + 2v (x) - 2v (y)]f j (x,y) = 0, (19) 

and Qj is a particular solution of (fTflj) . The only restriction on fj and Qj is that they should 
be selected in such a way that (x\Qj\y)* = (y\Qj\x), i.e., yield a Hermitian Qj. We can 
adopt, without loss of generality, a particular solution satisfying the Hermiticity condition 
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Qj(x,y) = Qj(y,x)*. In this case fj(x,y) = fj(y,x)* and as a result the general form of fj is 
independent of j. 2 

For the case that H = \p 2 , fj is given by d'Alemert's formula [To] : 

fj(x,y) = ^j(x -y) + xj(x + y), (20) 

where <Pj,Xj '■ —>■ C are arbitrary twice differentiate functions. 3 Because, according to (jl2j) . 
R 2 = 0, we can set Q 2 (x, y) = and use (fTHJ) and pUj) . to obtain 

(x\Q 2 \y) = ip 2 {x-y)+ X 2{x + y). (21) 

We can easily obtain the explicit form of the operator Q 2 using (f2*T]) . The result is Q 2 = 
Vzij*) + X2{p) where a tilde stands for the Fourier transform: F{j>) := e~ ipx F(x)dx, and 
V is the parity operator. The same reasoning shows that, whenever H = |p 2 , 

Q j = Q j +F j +K j V, (22) 

where j G Z + is arbitrary, Qj is a particular solution of the operator equation 0, and Fj = 
<fj{p) and Kj = Xj{p) are an y functions of p that together with Qj make Qj Hermitian. 4 In 
particular, if we select Qj to be Hermitian, Fj + Kj V will be Hermitian. 

It turns out that the condition that Fj + Kj V be Hermitian, which can be stated as 

VK]-KjV = Fj-Fl (23) 

is equivalent to the requirement that Fj and Kj be respectively Hermitian and "P-pseudo- 
Hermitian i.e., 

F] = F J} K] = VKjV' 1 = VKjV. (24) 

Therefore, as functions of p, Fj is real-valued while Kj is "PT-invariant. 5 One way of seeing 
this is to use the identities [p, V] = 2pV and \p, Kj] = to show that 

[p, V K) - KjV] = 2p{VK] - KjV) . (25) 

In view of (|2*3*jl and [p, Fj] = 0, the left-hand side of ()25|) vanishes identically whereas its 
right-hand side equals 2p(Fj — Fj). Hence, Fj — Fj = 0, which together with (^j) imply pijl. 

For the case that i?o involves a non- vanishing potential Vq, equations (JTB)) and ()19|) do not 
admit a closed form exact solution. In fact, one can prove that the factorization technique 
that is at the heart of the derivation of d'Alemert's formula fl2()jl does not apply to ()19j) . The 
spectral method that is also originally developed to solve wave (and heat) equations admit a 

2 We will however continue using the label j, because for different j we can adopt different homogeneous 
solutions. 

3 The Hermiticity condition: fj(x,y) — fj(y,x)* imposes the following restrictions on ipj andx^: 3?[<^j(— x)] = 
$t[tfij(x)], Ss[<pj(— x)] = —^s[ipj(x)] — 2rj, $t\Xj( x )] — r ji where 3? and 3 stand for real and imaginary parts of 
their argument respectively and rj is an arbitrary real constant. 

4 They must also be piece-wise analytic and possess (inverse) Fourier transforms. 

5 This means that for any real variable p c , Fj(j> c ) £ K and Kj(p c )* = Kj(—p c ). 
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generalization to both (fTBj) and t)19j). Yet it does not produce generic closed form expressions 
for the most general solutions of (fTBj) and (fTU|) . 

In the remainder of this paper we will only be concerned with the case v o = 0, i.e., H = \p 2 ■ 
In this case we can use the following prescription to compute Qj. 

(i) Fourier transform both sides of fTTj) over y to yield the ordinary differential equation: 

(d 2 x +p 2 )(x\Q J \p) = -2WM, (26) 
where p is treated as a free real parameter. This equation admits solutions of the form 

PX 

-2 / dx'Gpix^^ix'lRjlp) (27) 

J x 

where G p (x, x') = sin[p(x — x')]/p is the Green's function for d 2 + p 2 ; 

(ii) Evaluate the integral in (}2"Tj) by making a simple choice for x and identify the resulting 
particular solution of (f2*U|) with (x\Qj\p). 

(Hi) Take the inverse Fourier transform of (x\Qj\p) over p and denote the result by (x\Qj\y). 
If (x\Qj\y)* = (y\Qj\x) proceed to (iv). If not, add the homogeneous solutions u±(p) e ±ipx 
to (x\Qj\p), redefine (x\Qj\y) to be the inverse Fourier transform of (x\Qj\p)+u + (p) e ipx + 
u_(p) e~ ipx over p, fix u± in such a way that (x\Qj\y)* = (y\Qj\x), and relabel (x\Qj\p) + 
u + (p) e ipx + u^(p) e~ ipx by (x\Qj\p). 6 

(iv) Let Qj(x,p) := e~ txp {x\ Qj\p) and order the terms in Qj(x,p) in such a way that 
all x's are placed to the left of p's. This yields Qj, if x and p are identified with the 
corresponding operators. 

This prescription produces an expression for the operator Qj that is not manifestly Hermitian. 
To obtain a manifestly Hermitian expression, one must use the standard commutation relations 
to make the necessary reordering of the terms in Qj. A manifestly Hermitian expression that sig- 
nificantly simplifies the comparison of different operators (or different expressions for the same 
operator) is the one having the following symmetric form: go(x)+h (p)+Y^kLi Sk {gk(x), hk(p)}, 
where for all k e {0, 1,2, ■ ■ •}, and are real-valued nonconstant functions, 6 M, and 
{-, •} is the anticommutator, [3[7j. 

3 Metric Operators for v = iex 3 

The imaginary cubic potential, v(x) = iex 3 , has a real, positive, discrete spectrum [T6 | IT7 | IT8]. 
Therefore, according to the general results of pQ, it admits a nonempty set of positive-definite 
metric operators r/ + . Because this potential is purely imaginary, it is an ideal toy model to 
apply the general method developed in the preceding section. 

6 One can show that there is always a pair of functions u± that satisfy (x\Qj\y)* — (y\Qj\x). 
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The calculation of Qj for the potential v = iex 3 has been considered in ^3] where the 
authors obtain the operators Qj that yield the metric operator r/ + associated with the CVT- 
inner product. They achieve this in two steps: (1) They use the restriction imposed by the 
choice of the CPT-inner product to infer that Qj must be even in x and odd in p and to obtain 
a more restricted set of operator equations for Qj (This amounts to setting Qj = for even j.) 
(2) They adopt an appropriate ansatz for the form of Qj and determine its free coefficients by 
enforcing the operator equations for Qj. 

In this section, we will use our analytic method to obtain the most general Qj for the 
potential v = iex 3 . Our method is more systematic, technically more convenient, and more 
general, for it does not involve making any kind of (a priori) assumptions about the general 
structure or symmetries of Qj. In particular, it does not rely on making an ansatz for Qj. 
Indeed, we will see that the ansatz used in [T3j misses a large class of CPT-inner products even 
in the first order of perturbation theory. 

As we have argued in the preceding section for purely imaginary potentials such as v — iex 3 , 
the operators Qj have the general form (|22jl. In the following, we will first obtain the structure 
of the functions Fj and Kj appearing in (|2*2*|) and then present our calculation of Qj and Qj for 
J = 1,2,3. 

3.1 Calculation of Fj and Kj 

By definition [H EIJ El IE] , the physical observables O are linear operators that are Hermitian 
with respect to the inner product (-,-) v+ , i.e., they are ?7 + -pseudo-Hermitian operators [2]. 
We can construct them in terms of linear operators o that are Hermitian with respect to the 
reference inner product (-|-) according to O = or)+ 2 . For example, for o = x and o = p, 
we find the 77 + -pseudo-Hermitian position (X) and momentum (P) operators, 

X := V~ + 1/2 xvT = e Q/2 xe~^ = x + J2 [ -^r^QU (28) 

k=l 

P := V+ 1 * P tf = eWpe^ = p + £L±\p,Q) h , (29) 

k=l 

respectively, [U E] . These equations show how the operators Qj and in particular Fj and Kj 
enter the expression for the basic physical observables of the theory. 

Next, let I G K + and scale x according to x —>■ x/L The requirements that the canonical 
commutation relation is not affected by this scaling and that the Hamiltonian H = ^p 2 + iex 3 
changes by a total scaling (so that the scaling x — > x/i has the non-physical effect of changing 
units, as it should), we find that p — > £p and e — > £ 5 e. 

In light of (0), it is not difficult to see that the right-hand side of (J28|) involves terms of 
the form e^[x, Fj(p)] = ie^Fj(p) and e J [x, Kj(p)]V = ie 1 Kj(p)V, where a prime means the first 
derivative. Therefore, if we demand X to have the same scaling properties as x, i.e., X — > X/£, 
we find that x — > xjl implies 

Fj (p) - Fj (Ip) = i^Fj (p) , Kj (p) - Kj (lp) = i^Kj (p) . (30) 
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In particular, we have Fj(p) = p~^Xj(p) and Kj{p) = p~ 5 Wj(p) where Xj(£p) = Xj(p) and 
Oj(£p) = Oj(p) for all I G M + . Furthermore, the assumption that the quantum theory under study 
does not involve a preassigned hidden length scale implies that Xj and 9j must be constants. 7 
This together with the requirement that Fj be real-valued and Kj be PT-invariant then imply 
that Xj G K. and 9j = V Kj for some Kj G R. Therefore, for the imaginary cubic potential, we 
have for all j G Z + : 

FM = *)<p) = J*. (31) 

where Xj, Kj G M are arbitrary. 8 
3.2 Calculation of Qj 

For the potential v = iex 3 , we have Hq = \p 2 ■, Hi = ix 3 , and R\ = —2ix 3 . Hence, the most 
general form of Qi is obtained by setting j = 1 in (|22j). where (x|<2i|?/) is a particular solution 
of the wave equation 

[-dl + ^]<z|Qi|y) = -Azx 3 5(x - y), (32) 

and 5(x) denotes the Dirac delta function. To obtain Q\ we use the prescription described in 
Section 2. Taking the Fourier transform of ()32|) and solving (|26|) for a particular solution, we 
find 

(x\ Qi b) = e ip:c f x 4 ^ + 2i a; 3 ^ - 3 a; 2 ^ - 3* x J . (33) 



Hence, 



2V27T VP P P P' 

1 / ,1 _. , 1 . ,1 _. 1 



Qi = ~(x 4 - + 2ix 3 — -3x 2 ^-3ix^r) . (34) 
2 \ p p A p 6 £T / 

We can evaluate the inverse Fourier transform of (J33)) to obtain 

(x\Qi\y) = -xy(x 2 + y 2 )sign(x -y), (35) 

where sign(x) := x/\x\ for x^O and sign(O) := 0. 

Eq. (|33j) shows that indeed Q\ is a Hermitian operator. With the help of the identities 

[x, f(p)} = i/'(p), [x 3 , f(p)\ = | {x\ f'(p)} + ~ /< 3 )(p), (36) 

where stands for the n-th derivative of / and /' := /W, we can express Qi in the following 
manifestly Hermitian form 

2i = l{^-} + |{* 2 ,4} + 4r ( 37 ) 
4 p 4 p d p a 

A useful check on the validity of our calculation of Q\ is to verify that ()35|) satisfies (|32|). 
We have checked this by direct substitution of in the left-hand side of 



7 This argument is effectively equivalent to and provides a conceptual interpretation for what the authors' of 

refer to as "dimensional consistency" . 
8 We can also apply the assumption of lack of a hidden length scale in the underlying classical theory to 
arrive at fl5T]). 
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Substituting (|37jl in (j22j) with j = 1 and using (JBTj) we find the most general form of Qi, 
namely 

^ 1 r J h 3 , n 1 -. Ai iKl _ , . 

Ql = T -} + 7 -} + ^ + "If P ' 38 

4 p 4 p d p b p b 

where Ai := Ai + 3 and K\ are arbitrary real constants. 

Now, we are in a position to compare our result for Qi with that obtained by Bender, Brody, 
and Jones (BBJ) ^3], namely 

QBBJ = / 4i + 4x 3l ^ + g x 2l x 2 + 4x 1 ^3 + 1 A + « ,g g N 

32 VP p p p p / p 

where aGlis arbitrary. We can use the identities, 

x 3 -x + x-x 3 = {x\ -} + 3{x 2 , — } + — , x 2 - x 2 = - {x\ -} + 2{x 2 , — } + (40) 
p p p p A p° pip p A p° 

to express (|3^jl in the form 

Qf BJ = j{^-} + |{x 2 ,4} + (« + x)4- ( 41 ) 

4 p 4 p J 4 p & 

Clearly, Q^ BJ defines a one-parameter subfamily of the operators of the form (|38jl . It corre- 
sponds to the choice Ai = a + and Ki = 0. 

It is interesting to see that (|3*Hj) does also define a CPT-inner product, for it satisfies all the 
conditions stated in [TJ]. In particular, if we use the number operator N = |(x 2 + p 2 — 1) to 
express V in the form 9 

V = (-If = e i7rJV = e^ 2+x2 ^ = -ie^ { v 2+x2 \ (42) 

we see that the last term in (J38|) is also even in x and odd in p. The analysis of ^3] seems to 
have missed this term, because it relies on the choice of a particular ansatz for Q\. 

Our derivation of Qi shows that, under the assumption that the quantum system defined by 
the Hamiltonian H = y + ie x 3 does not have a hidden length scale, up to (and including) the 
first order terms in the perturbation parameter e, the CPT-inner products are the most general 
inner products that restore the unitarity of the dynamics. This result does not, however, extend 
to higher orders in perturbation theory. It already fails in the second order, because in general 

Q 2 = F 2 (p) + K 2 (p)V = 4 (A 2 - K2V), (43) 

whereas for a CPT-inner product Q2 = 0. This is a manifestation of the fact that the CPT-inner 
products form a proper subset of the set of all allowed inner products. 

Next, we wish to compute Q3 using our method. We postpone the details of this calculation 
to Appendix B. Here we outline its general strategy. 

9 Recall that the eigenfunctions (x\n) of the harmonic oscillator Hamiltonian \{x 2 +p 2 — 1) are even (respec- 
tively odd) functions of x for even (respectively odd) values of the spectral label n. Hence V = (— 1) . 
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The first step in the calculation of Q3 is the determination of R3 which in view of f!13|) is 
given by 

R3 = -\[[x\QilQi]. (44) 

We recall that according to the prescription explained in Section 2 we need to compute (x\Rz\p). 
Therefore, instead of evaluating the double commutator in (jUJ), we compute 

(x\R 3 \y) = --(x\[[ix 3 ,Q 1 ],Q l ]\y) = -- / dz(x 3 + y 3 - 2z 3 )(x\Q 1 \z)(z\Q 1 \y), (45) 
6 6 

and take its Fourier transform over y. 10 Next, we insert Q\ = Q\ + \p~ 5 + inp~ 5 V, which is 
equivalent to (jSHJ), in Eq. (f4*H|) and use to perform the integral over z. Doing the necessary 
calculations, as outlined in Appendix B, we then find 

R3 = 5*0,0 + AiS^o + ftiSo.i + AiftiSi,! + A^S^o + 5*0,2, (46) 

where the operators S^ u are defined by Eqs. (J%U|) - (J%o"|) below. 

In view of f!46|) . we can construct a particular solution of (J7|) of the form 

Q 3 = T 0i0 + AxT^o + KiT 0>1 + \iK X T x>1 + \\T 2fl + ^T , 2 , (47) 

where T^ u solves the operator equation 

[p 2 ,T^]=2S^. (48) 

Hence, (x\T^ tl ,\y) is a particular solution of the wave equation 

{-% + %)(x\T lv ,\y)=2(x\S lv ,\y). (49) 

In Appendix B we construct such solutions and check that the corresponding operators, T^ v , 
are Hermitian. This shows that the operator Q3, as given by (|47jl . is Hermitian. Consequently 
Eqs. (jSU) hold for j = 3. 

Using (|2*2*jl . (|3*T|) . (JUJ), and the explicit form of the operators given in Appendix B, we 
obtain after a lengthy calculation 

5 7 / 4 _ 



q 3 = E ^ {** ^} +^-i(J>* ^> + ^ ) ^ ( 5 °) 



where 



^01 := Cool + Ai C101 + A^C201 + «lC021j ^02 := C002 + Ai C102 + K1C012 

<^03 := C003 + Ai C103, G?04 = 0)04, C?05 = 0)05 ^11 := ^i(con + AiCm), 

di2 ■= Ki(coi2 + A1C112), dia := K1C013, du := K1C014, 

A3 := A3 + 2(cqoo + A1C100 + A^C200 + ^lCo2o), ^3 := K 3 + 2/ti(Coio + AiCno), 
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The advantage of this approach is that the above calculation can be done using Mathematica. 
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c^e are real constants given in Tabled Ai, K\ are the free real parameters determining Qi, and 
A3, K3 (or equivalently A3, k 3 ) are another pair of free real parameters. 

Note that the freedom in the choice of parameters a and /3 entering the expression obtained 
in pi] for Q3 corresponds to the freedom in the choice of Ai and A3 (or A3). Furthermore, it is 
not difficult to see that Q3 constructed above fulfils all the conditions imposed by the choice of 
CPT-inner product. The fact that the approach of ^3] do not reveal the presence of the terms 
proportional to K\ and £3 stems from the particular choice for the ansatz used in [14] . 

We conclude this section by mentioning that we can similarly apply our method to obtain 
the general form of Qj for j > 3. But as expected the algebra becomes quite involved. 

4 Physical Implications of Changing the Metric Opera- 
tor in Quantum and Classical Treatments 

In order to understand how the freedom in the choice of the metric operator affects the physical 
content of the theory, we will use the most general metric operator to perform a perturbative 
calculation of the pseudo-Hermitian position (|2*8*|) and momentum (|29jl operators, the equivalent 
Hermitian Hamiltonian, 

h-^^Hr,- 112 = e- Q l 2 He Q '\ (51) 
and the classical Hamiltonian 11 , 

H c (x c , p c ) := hmh(x,p) , (52) 

ft- »0 x-+a; CJ p— »-p c 

for the system defined by the imaginary cubic potential. 
We can easily show, using ©, 1|2B|) . flZfl and 0, that 

X = x-i[a;,Q 1 ]e + 0(e 2 ), P = p - l - [p, Q x \ e + 0(e 2 ), (53) 

h = ^+ i -[x\Q 1 }e 2 + + l -ix 3 ,Q 2 }e 3 + 0(e% (54) 

where 0(e n ) stand for terms of order n and higher in powers of e. Therefore, we first use f)38jl 
and (|3HJ) to compute 

[x , Qi] = -\ Ux\ 1} + 9 {x 2 , 1} + ^ - 4* {x, 1} Vj , (55) 

ip,qi) = 4( 2{x3 ^ }+3{ v } -^ p )' (56) 

[x\Q x ] = _^('{x 6 ,l} + 22{x 4 ,l} + (510 + 10A 1 ){x 2 ,4}+ 
L J 4 \ p z p 4 p b 



8820 + 140A 



I-l^a* (57) 



p8 3 p5 

( 1 44 \ 1 
[x\Q 2 ] = -15zA 2 ^ 2 , — } + — J - K2 {a; 3 , — } V. (58) 

a; c and p c appearing in l|52|l are respectively the classical position and momentum observables. 
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Substituting these relations in and (p>4"j) . we have 

A' , +- I {,-._} + 9{x 2 i} + ^-4 Kl {x ,l}v\e + 0(e 2 ), (W) 




P = p+ -(2{x 6 ,-} + 3{x,-}-^V )e + 0(e 2 ), (60) 
4 \ p p 6 p 4 J 

h = ^ + ^f{x 6 ,l} + 22{x 4 ,i} + (510 + 10A 1 ){a: 2 ,^}+ 
2 lo \ p° 



8820 + 140Ai 4 r , 1 \ , 
p$ 3 p i 

\ (^\ 2 ({x\ ^} + ^) ~ ^} Vj + 0(e 4 ). (61) 

Equations (|59|) and (|60|) show how the free parameters Ai and k%, that determine the metric 
operator up to terms of order e, enter the definition of the basic observables of the theory. 

The quantum theory defined by the imaginary cubic potential can be described by the 
manifestly Hermitian Hamiltonian h within the framework of the standard quantum mechanics, 
[27H ITH] . However, as seen from Eq. (|6*Tj) this equivalent Hermitian description is sensitive to the 
choice of the metric. The metric-independence of the Hamiltonian H and the metric- dependence 
of the physical observables O in the non-Hermitian description of the quantum system are traded 
with the metric- dependence of the Hamiltonian h and the metric-independence of the physical 
observables o in its Hermitian description. 

Having calculated the Hermitian Hamiltonian h, we can determine the classical Hamiltonian 
(p)2*j) . This requires making the ^-dependence of the terms in (joT|) explicit. We do this by 
letting h — > £~ 2 m~ 1 h 2 h, x — > £~ 1 x, p — > ifr x p, where i is an arbitrary length scale and m is the 
mass, [Hj. Making this transformations, replacing the quantum observables with their classical 
counterparts, i.e., x — > x c and p — > p c , and taking the limit h — > 0, we find the following 
remarkably simple expression for the classical Hamiltonian (J52j) . 

H c = f + lme 2 ^ + 0(e*). (62) 
2m 8 p z c 



Clearly H c is an even and nonnegative function of x c and p c . This is an indication that it 
supports closed classical phase space orbits. As shown in Figure 1, this is actually the case. 

A more important observation is that the terms in (JBTj) that involve Ai, K\, X2 and K2 do 
not contributes to the classical Hamiltonian, because they involve positive integer powers of h 
that vanish in the classical limit h — > 0. It is not difficult to see that the same behavior holds 
in all orders of perturbation; the terms involving Xj and Kj, that characterize the freedom in 
the choice of the metric operator, do not contribute to the classical Hamiltonian. This means, 
at least for the system we consider, that the classical Hamiltonian is not sensitive to the choice 
of the metric operator. We expect this assertion to hold true generally. 
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-2-10 1 2 

x c 

Figure 1: Typical phase-space orbits for the classical Hamiltonian H c , for e = 0.1 and m — 1. 

5 Pseudo-Hermitian Quantization of the Free Particle 

If we set H = p 2 /2 and v o = Hi = in the constructions of Section 2, we find the most general 
metric operator t] + = e~® that renders the free particle Hamiltonian H = Hq Hermitian. It is 
not difficult to see that in this case Rj = for all j G Z + and that 

Q = F + KV, (63) 

for a pair of functions F and K of p. Again the fact that the system does not involve a particular 
length scale implies that F(p) =: A and K(p) =: k are real constants. Hence 

r]+ = e- x e- KV = e- A [cosh(K) - sinh(«)P]. (64) 

For this system the pseudo-Hermitian position (|2*5|l and momentum (|2~9|l operators take the 
form 

X = x e~ RV = e KV x, P = p e~ Kp = e Kp p. (65) 

Clearly the Hermitian hamiltonian h and the classical Hamiltonian H c coincide with H and 
pl/2, respectively. The pseudo-Hermitian quantization scheme defined by 

x c — > X p c — ^ P, Poission bracket — > — i commutator, (66) 

is the most general physically admissible canonical quantization of the free particle in one 
dimension that is compatible with the usual form of the free particle Hamiltonian. Note that 
for k^O the operators X and P are not Hermitian with respect to the usual L 2 -inner product. 
Yet they furnish a unitary irreducible representation of the Heisenberg-Weyl algebra [x,p] = i. 
The representation space is the Hilbert space 7^ p hy S = : ^ ~~ * < £\(4 , ,4>)ri + < oo} endowed 
with the inner product, 

/oo poo 
dx 4>(x)*ip(x) — e _A//2 sinh(K) / dx 0(x)*^(— x), 
-oo J — oo 

(67) 
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where ip, : R — ► C are arbitrary functions. 12 

As suggested by (J55J) and (|67j). the pseudo-Hermitian quantization of the free particle may 
be viewed as a smooth deformation of the standard (Hermitian) canonical quantization (which 
corresponds to taking k = 0). 

Next, we recall that 77+ defines a unitary operator mapping 7i p hy S onto L 2 (R), [201 El- 
Under this mapping X and P are transformed to x and p respectively and H is left invariant. 
Therefore, the pseudo-Hermitian quantum systems associated with all A, k G R are physically 
equivalent. In particular, there is no physical reason to choose k — 0. We should like to warn 
however that taking k^O has unusual yet completely consistent consequences. For example, 
the localized states are described by wave functions that in general have two delta function 
singularities! Up to an irrelevant phase factor, the delta-function normalized localized state 
vector centered at y G R is given by := f] + 1 ^ 2 \y) 1 jSHTij Therefore, according to ([fiijl . 

( x \^y)) = e A/2 [cosh(f ) 8{x -y) + sinh(f ) 5{x + y)\. 

Note however that this is not the physical position wave function for the localized state. The 
latter has the form {^ x \ ^ y ') v+ = ( x \y) — — v)- 

Finally, we would like to point out that according to (|65|). X 2 = x 2 , XP = xp, P 2 = p 2 . 
Therefore if we perform the pseudo-Hermitian canonical quantization (jfifij) on any classical 
Hamiltonian H c that can be constructed out of x 2 , x c p c and p 2 , we obtain the same expression 
for the quantum Hamiltonian as the one obtained using the usual canonical quantization. 
Conversely, every quantum Hamiltonian H that can be expressed in terms of x 2 , xp + px, and 
p 2 admits a metric operator of the form Ij64j) . 



6 Higher Dimensions and Scaling-Invariance 

The situation is considerably more complicated in n-dimensional configuration spaces with 
n > 1. In this case the operator equations (J7|) correspond to certain ultrahyperbolic equations in 
(n + n)-dimensions. These equations have peculiar properties, and except for very special cases 
their standard Cauchy problem is not well-posed {21]. For the free particle in n-dimensions, Q 
satisfies 

(-V^ + V 2 )(f|Q|y}=0, (68) 

where x = (xi,x 2 ,--- ,x n ) and V 2 := Y17=i^x - ^ i s no ^ difficult to see that does not 
provide the most general solution of ([61^) . 13 It only provides a special class of solutions. An 
important distinction between this class for n > 1 and the corresponding general solution for 
n — 1 is that for n > 1 the requirement that the theory must not depend on an arbitrary length 
scale does not restrict F and K to constants. This is simply because one can now use the ratios 

12 Bccause the metric operator (|64|) is a bounded operator acting in L 2 (R), 7i p hy S and L 2 (R) are identical as 
vector spaces. What distinguishes them from one another is their inner products. 

13 For example, for n = 2, the angular momentum operator, Q — X1P2 ~ which is clearly not of the form 
(|63|l . provides a solution of (|68(l . 
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of the components pi of the momentum operator p to construct scaling-independent operators, 
i.e., require that F and K are functions of the homogeneous p-coordinates. For example for 
n = 2, we can take F = F(pi/p 2 ) and K = K(pi/p 2 ). 

A similar situation occurs in one-dimension provided that the system involves an e-independent 
length scale. A simple example is the scattering potential, 

{— iesign(x) for \x\ < |f 
(69) 
for \x\ > § , 

which includes the length scale L. This potential has a real spectrum [Jj. It is also imaginary. 
Therefore, we can apply our method to obtain the general form of Qj. In j7] we have used a 
considerably more difficult spectral method to obtain the following particular solution of (fT7j) . 

Qi(x,y) = {x\Qx\y) = ^(\x + y + L\ + \x + y - L\ - 2\x + y\ - 2L) sign(a; - y). (70) 



According to the results of Section 2, the general form of Qi is given by (|22j) with j — 1. 
However, unlike the case of imaginary cubic potential, the operators Fj and Kj appearing in 
fl22(l can now be complicated functions of Lp. Specifically, given the fact that for the system 
defined by the potential (|69J1 . the scaling transformation x — > x/t induces the transformation 
e — > £ 2 e, it is not difficult to infer that Fj and Kj have the form 

J p 2i ' p 2j ' 

where Aj and Kj are respectively real- valued and PT-invariant functions. In particular, we have 

Qi = Qi+p- 2 [\i(Lp) + Kl (Lp)V], (71) 
where Q\ is given by (|70|) and 

Q 2 = p- 4 [A 2 (Lp) + K2 (Lj9)P]. (72) 

The derivation of (|72jl using the spectral approach pursued in [7j is an extremely difficult (and 
open) problem. 



7 Concluding Remarks 

In this paper we offer a method to compute the most general metric operator for a given quasi- 
Hermitian standard Hamiltonian in one-dimension. This method can be conveniently used 
for imaginary potentials. It allows a better understanding of the issue of non-uniqueness of 
the metric operator and the positive-definite inner products that render the quantum theory 
unitary. 

Applying this method to the imaginary cubic potential we discovered a new class of CVT- 
inner products and established the existence of admissible non-CPT-inner products. Further- 
more, we showed that the freedom in the choice of the metric operator affected the basic physical 
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observables of the theory in its non- Hermit ian representation and the Hamiltonian in its Her- 
mitian representation. An interesting consequence of our investigation is that the underlying 
classical system is not sensitive to the choice of the inner product. 

We also used our method to obtain a complete characterization of pseudo-Hermitian quan- 
tization of the free particle Hamiltonian in one-dimension. 

We wish to emphasize that the choice of the metric operator can be related to that of 
an irreducible set of observables as explained in [22] • Specifically, the requirement that the 
members of an irreducible set of operators 14 be Hermitian fixes the metric operator (up to an 
irrelevant constant coefficient) [221 [2B]- The application of this program for "PT-symmetric 
Hamiltonians having a real spectrum involves selecting sufficiently many operators (all having 
real spectra) such that together with the Hamitonian they form an irreducible set. These 
operators must however be compatible in the sense that there must exist a positive-definite 
metric operator r] + such that all of them be ?? + -pseudo- Hermitian. In general the determination 
of the compatibility of the operators that are to be chosen is a difficult task. In particular it 
requires the knowledge of the most general positive-definite metric operator that renders the 
Hamiltonian pseudo-Hermitian. Therefore, the results we have reported in this paper play a 
central role in employing the method of [22] i n PT-symmetric quantum mechanics. 

End Note: For the imaginary cubic potential considered here it is sometimes argued that 
(x\r)\y) is not an analytic function of e and its series expansion studied in ^3] and here is 
meaningless, because e is not dimensionless and scaling x and y one can make the numerical 
value of e as large or small as one wishes. This argument is inconclusive, for it certainly applies 
to e ex5 , which is an analytic function of e for all x G M, and to (1 + ex 5 )" 1 which is analytic 
in e for \x\ < |e| 1//5 . As these examples suggest, in general there is a region in the x-y plane in 
which (x\i]\y) admits a power series expansion in e. The issue of the summability of this series 
is an extremely difficult open problem. The fact that e is not dimensionless implies that an 
approximate evaluation of rj that is based on a truncation of the power series for (x\r)\y) will 
be valid only within a region T> e of the x-y plane whose size depends on e. Obviously, because 
the contribution of every order in e to rj is dimensionless, scaling e would scale the size of T> e . 



14 A set of operators acting in a Hilbert space TL is called irreducible if there is no proper subspace of TL that 
is left invariant under the action of all members of this set. 
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Appendix A: Derivation of the Operator Equations Sat- 
isfied by Qj 

In order to derive (J7J) we first substitute (jHJ) and (jlj) in (j2J) and use flBJ) and the fact that 
H{ = —Hi to obtain 

etf 1 = -i ({H ,Ql + leH u Ql) (73) 

where 

[A, Ql := e^Ae* - A = - [A, Q] fc . (74) 

fe=i 

Next, we substitute the right-hand side of Eq. (|75|) for the eH\ in the term [eifi,<5]i that 
appears in this equation and do the same in the right-hand side of the resulting equation. 
Repeating this procedure t — 1 times yields 

^i = 2EE *«*l H °> Ql + [eH u Ql, (75) 

m=l fc=m 

where i G Z + is arbitrary and for all k, m G Z + 



m 



n=l 



fc!2 m 1 \ n \ n I n\(m — n)\ 



In view of (j3J), (J74*j) . (|75)l. and (|75jl . it is easy to see that 



£ I 



tH ^ = \HH «™*[#o, Q]* + (77) 

m=l k=m 

where 0(e e ) stands for terms of order i or higher in powers of e. We can employ the identity 



Em=i ELm = £k=i £m=i to express dZZD in the form 



eF 1 = -X>[#Q,QL+^ +1 ), (78) 



k=l 



where g fc := E m =i 9mk- 

Next, we substitute (j3J) in(J7SJ) and require that (|7Hj) be satisfied at each order j < £ of the 
perturbation. For j = 1, this implies 

^ — [tfo.Qi], (79) 

where we have used q± = — 1. For j > 2, we solve for e*[Ho, Qj] in terms of the remaining terms 
of order e J in (fTSj) . This together with (J79J) yield (JJJ) . 
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Appendix B: Calculation of Q3 

In view of flU} , the operator S^ u , that determine R3 according to (HHJ), are given by 



\X 



\ x 



U' 



U' 



[ x 



\Sq,o 

1^2,0 
1^0,2 



\y) 
\y) 
\y) 
\y) 
\y) 
\y) 



% 

i 

■y roc 

6 
1 

6 



-00 

00 



dz{x 3 + y 3 - 2z 3 ) (xlQ^z) (z\Q x \y), (80) 

dz(x 3 + y 3 - 2z 3 )[(x\Q 1 \z)(z\p- 5 \y) + {x\p- s \z) (z\Qx\y)], (81) 

dz{x 3 + y 3 - 2z 3 ){(x\Q 1 \z)(z\p~ 5 \ -y) + (x\p- 5 \ - z)(z\Q 1 \y)}, (82) 

dz(x 3 + y 3 - 2z 3 )[(x\p~ 5 \z)(z\p~ 5 \ - y) + (x|p~ 5 | - z) (z\p~~ 5 \y)}, (83) 

-00 

dz{x 3 + y 3 -2z 3 )(x\p- 5 \z)(z\p- 5 \y), (84) 

^(x 3 + y 3 -22 3 )(x|p~ 5 | -z){z\p^\ -y) (85) 



-00 
00 



-00 
00 



We can use (J33|) and the identities 



sign (x J 



1 ^00 ixk 
- dk 



ITT 



k: 



(x\ — \y) = — dk 



P 



2tt 



k 5 



— (x -y) 4 sign(x -y), (86) 



to evaluate the integrals appearing in (|8U|) - (|85|) by turning them into a pair of Fourier trans- 
forms that we can easily perform using Mathematica. The result is 



{x\So fi \y) 

(x\s ,i\y) 

(x\Si,o\y) 
(x\Si,i\y) 

(x\S ,2\y) 

(x\S 2 ,o\y) 



26880 
1 



319334400 



[l5(x n y - x y 11 ) + 7(x 9 y 3 - x 3 y 9 ) - 48(x 8 y 4 - x 4 y 8 )] sign(x - y), (87) 
(x + yf [3115(x 7 - y 7 ) - 3818(x 6 y - xy 6 )+ 



6386688 
1 

47900160 

i 

95800320 

i 

6386688 



3120(x 5 y 2 - x 2 y 5 ) - 12123(xV - x 3 y 4 )] sign(x + y), 
(x - yf [623(x 5 + y 5 ) + 1970(x 4 y + xy 4 )+ 

3743(x 3 y 2 + x 2 y 3 )} sign(x-y), 
(x + yf [7{x 3 - y 3 ) - I5(x 2 y - xy 2 )} sign(x + y), 

(x - yf [29(x 3 + y 3 ) + 15(x 2 y + xy 2 )} sign(x - y), 

(x - yf 1 (x + y) sign(x - y). 



(90) 
(91) 
(92) 



These equations show that, as expected, S^ u are anti-Hermitian. 

Next, we use (J%7|) - ((02]) and the the prescription given in Section 2 to transform the wave 
equations (J4T?|) into the ordinary differential equation (d 2 +p 2 )(x\T^ v \p) = — 2(x\S fl ^\p), obtain 
particular solutions of this equation, and determine the corresponding operators T^ u that solve 
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. This yields 
To,o 



10 



1,0 



-0,2 



£=1 ^ £=1 



aoi£ £ ^15^1 T 3 , 



1 



P 



V, 



£=1 ^ 



2 1 



(93) 
(94) 
(95) 



where are real and positive constants given in Table El 

Next, we wish to determine if T^ v and consequently the operator Q 3 given by (}4Tj) are 
Hermitian. We can easily answer this question by inspecting the matrix elements: 



x 



x 



X 



X 



X 



Tq,i\p) 

Ti,i\y) 

T 2,o\y) 



i(x — 



[X 



l[X 



X 



y) 4 £&oo* (xV^ + ^-V) sign(x-y), 
t=\ 
4 

+ yf bou (xV- 1 + x 8 ~ e y e ) sign(x + y), 
3 

- yf b m (a?V~' + x 6 - e y £ ) sign(a; - y), 
2 

+ y) w b iu (*V~' + x4 ~V) sign(x + y), 



£=1 
2 



i(x 



y) 10 6 °2£ OV t + x A V) sign(a; - y), 



i b 20 i (x ~ y) xy sign(x - y), 



(96) 

(97) 

(98) 

(99) 

(100) 
(101) 



where b^t are the real coefficients listed in Table El Eqs. (jHSj) - (jlOip show that indeed T^ v 
and consequently Q 3 are Hermitian operators. Hence in view of (j22J), (EH) and (j4*7|) . Q 3 has the 
following general form 

Q 3 = T ,o + AiTi )0 + «iT 0i i + AiKiTi,! + A 2 T 2 , + k\T^ 2 + A 3 p~ 15 - iK 3 p~ 15 V, (102) 

where Ai and k% are the real free parameters that fix Q\ and A3 and k 3 are a pair of arbitrary 
real constant coefficients. 

In order to derive a manifestly Hermitian expression for Q 3 , we first derive, after a lengthy 
calculation, the following manifestly Hermitian form of T^ v . 



Yll=0 c nvt {x 2£ , 15-21 } 



for v 7^ 1 



T, 



(103) 



- i Y.i=Q c ^u{x 1 \^2i]V for z/=l, 
where c^g are real constants given in Tabled Substituting ()103|) in ()102|) we finally find ()50|). 
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i 





1 


2 


3 


4 


5 


C 0(M 


141274966833 
32 


3830434839 
64 


23858793 
64 


43479 
32 


267 
64 


1 

80 


cqu 


24081603 
20 


328947 
40 


16327 
80 


35 
48 


1 

480 




Cioe 


54563145 
16 


1430535 
16 


8695 
16 


5 
3 








cue 


1547 
4 


61 

4 


1 

12 










CQ21 


-357 


9 


1 

12 











C2oe 


-2275 
4 


-25 
4 















Table 1: Numerical values of the coefficients c^e 
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1 


2 


3 


4 


5 


6 


7 


8 


9 


10 




2745171 
32 


2745171 
32 


677457 
16 


439857 
32 


52029 
16 


9375 
16 


651 

8 


273 
32 


5 
8 


1 

40 




70317 
5 


23592 
5 


20207 
20 


794 
5 


777 

40 


217 
120 


7 

60 


1 
240 








1110915 

8 


363315 
8 


36355 
4 


9305 
8 


90 


10 

3 












351 
2 


71 
2 


11 

3 


1 
6 
















388 


48 


11 

3 


1 
6 














o-2oe 


325 
2 


25 
2 



















Table 2: Coefficients a^e of the operators 7) 
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1 


2 


3 


4 


5 


booe 


79 


79 


533 


947 


53 


11468800 


2867200 


8601600 


8601600 


983040 


boie 


601 


4757 


5443 


937 




532224000 


1596672000 


798336000 


266112000 


bm 


211 


533 


9127 






18923520 


63866880 


510935040 


bin 


1 


1 








70963200 


383201280 


b()2£ 


13 


15 








479001600 


958003200 


b2oe 


1 

76640256 











Table 3: Coefficients b^t of the matrix elements (x\T^ v \y) 
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